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I. INTRODUCTION 

Collective transport effects in ordered many-fermion 
and many-boson systems include some of the most dra- 
matic and profound phenomena that occur in condensed 
matter physics. For example, the (practically) dissipa- 
tionless transport of electrical charge by Cooper pairsi in 
superconductors and superfluidity in 4 He and 3 He have 
been important topics through most of the field's his- 
tory. New instances of this general class of phenomenon 
continue to arise and create interest. One recent case is 
collective charge transport in double-layer quantum Hall 
systems, in which spontaneous inter-layer phase coher- 
ence leads to a strongly enhanced zero-bias tunnel cur- 
rent from one layer to the otheri2iiM£ Closely related 
issues connected with the possibility of superfluidity due 
to excitonic Bose condensation in electron-hole double- 
layer systems^ also continue to attract attention^ and 
inspire experimental activity. 

Two of us~ have recently proposed the possibility of 
realizing nearly dissipationless collective spin currents in 
easy-plane thin- film ferromagnets. Although the pro- 
posed phenomenon has some formal relationship to su- 
perconductivity and superfluidity, it depends in part on 
a symmetry (magnetization orientation invariance within 
an easy plane) that can only be realized approximately. 
More importantly, unlike the case of superconductors 
which can easily be biased by a current source, strate- 
gies for driving a real thin-film ferromagnet into the 
metastable spin-current state present experimental and 
materials challenges that have not yet been overcome. 
For ferromagnets, the meta-stable states that have a non- 
zero spin current are spiral states, as sketched in Fig. ^ 

In Ref. we demonstrated that these states are 
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FIG. 1: State with spiral magnetic order in the x-y plane 
characterized by the wave vector Q. The arrows represent 
the local spin density of the itinerant-electron system. 



metastable only when the system has a non-zero mag- 
netic anisotropy that is accurately uniaxial with an easy 
plane, and proposed an experimental strategy to gener- 
ate and detect macroscopic dissipationless spin currents 
in thin-film ferromagnets. Other related recent work has 
addressed persistent spin currents in rings that experi- 
ence an inhomogeneous magnetic field >i2iiiii£ In this pa- 
per we consider the case in which the ground state of 
the system has spiral order, the case of helimagnets. As 
mentioned already in Ref. |9|, the persistent spin current 
of a spiral state is proportional to the derivative of its 
energy with respect to the state's winding wavevector 
Q, a quantity which necessarily vanishes in the ground 
state. Like ferromagnets, helimagnets do not normally 
carry a spin current in their ground state. However, as 
we demonstrate in this paper, in the helimagnet case a 
persistent ground state spin current is induced by an ex- 
ternal magnetic field. 

Throughout this paper we consider states in which the 
spin density (s(r)} lies in the x-y plane. The stabil- 
ity of these states will, in general, require some easy- 
plane anisotropy that could have magnetostatic or mag- 
netocrystalline origin. The presence of such a magnetic 
anisotropy energy term is implicitly assumed through- 
out the paper. In the first part of the paper, we set the 
applied magnetic field to zero and consider states with 
simple spiral order, illustrated schematically in Fig. ^ 
Our objective in this section is to demonstrate by an ex- 
plicit microscopic calculation the property that the spiral 
ground states of helimagnets, unlike the spiral metastable 
excited states of ferromagnets, do not carry persistent 
spin currents. The order parameter of a spiral state is 

(s(r)) Q = s Q [cos(Q • r)x T- sin(Q ■ r)y] . (1.1) 

and it is characterized by a wave vector Q, by the mag- 
nitude sq of the order parameter, and by a chirality. For 
a single-band lattice model < sq < 0.5n c , where n e is 
the electron density per atomic site. As Q is varied, the 
change of the order parameter's spatial dependence will 
cause a change in the magnetic condensation energy and 
in the magnitude of the order parameter. We have ar- 
gued previously^ that these systems carry a dissipation- 
less spin current, i.e., a current with equal magnitude 
and opposite directions for up and down spins, that is 
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related to the dependence of the total energy density e 
on Q by 9 , 
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(1.2) 



Note that the spin quantization axis is along the z- 
direction here. All spin currents discussed in this paper 
are polarized along the direction perpendicular to the 
plane defined by the spiral order parameter; we therefore 
do not explicitly indicate the tensor character of spin cur- 
rents. Eq. Ql. 2D implies that it is not the presence of spiral 
magnetic configurations, but rather the change of total 
energy with varying spiral wave vector Q which leads to 
spin currents 13 . If the energy has a local minimum at a 
finite wave vector Q* (as in a helimagnet) then this state 
will, unlike spiral states in a system with a ferromagnetic 
ground state, not support a persistent spin current. We 
examine this property from a microscopic point of view 
in Section II. In Section III we discuss persistent spin cur- 
rents using the picture of "Cooper pairs" familiar from 
the Bardeen-Cooper-Schrieffer (BCS) mean-field theory 
of superconductivity, which provides an alternative way 
of understanding when and why spin supercurrents ap- 
pear for a given magnetization configuration. 

With this background, we turn our attention in Sec- 
tion IV to a new strategy for generating persistent spin 
currents. We point out that persistent spin currents can 
be generated in systems with spiral magnetic order sim- 
ply by applying a magnetic field oriented in the plane of 
spiral order, the x-y plane in our notation. The external 
magnetic field creates a competition between spiral order 



with wave vector Q* and homogeneous magnetization 
along the field direction, leading to magnetic order de- 
scribed by a soliton lattice with wave vector |Q| < |Q*|. 
We demonstrate that soliton-lattice states carry persis- 
tent spin currents, whose amplitude is controlled by the 
strength of the applied magnetic field. Finally, we con- 
clude in Section V with some brief speculations on possi- 
ble applications in spintronics for these field-induced spin 
currents. 



II. MICROSCOPIC MODEL FOR PERSISTENT 
SPIN CURRENTS IN HELIMAGNETS 

A. Hartree-Fock theory for spiral magnetic order 

In this section we try to shed additional light on the 
relationship between persistent spin currents and the de- 
pendence of energy on the spiral wavevector by carrying 
out explicit calculations for a simple toy model of an 
itinerant-electron system with magnetic order. Similar 
spiral state models have been presented previouslyi^i^; 
the following brief description is included for complete- 
ness and intended to establish notation for the following 
discussion. 

We consider a system of fermions on a lattice with 
single-particle band energy ek and delta-function repul- 
sive particle-particle interaction U8(vi — Yj), which we 
treat in a mean-field approximation. The unrestricted 
Hartree-Fock mean-field Hamiltonian for the spiral or- 
dered state with wave vector Q is 
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where h = Us Q and s Q = (l/V0Ek( c k+Q/2,T c k-Q/2,l)- 
The off-diagonal terms in Eq. H2.lt couple electronic 
states with opposite spin and Bloch wave-vector differ- 
ence Q. (The wavevectors above are understood to be 
reduced to the first Brillouin-zone of the lattice.) To di- 
agonalize the Hamiltonian we employ the transformation 
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with tan 20k = 2/i/[ek+q/2 — e k-Q/2] an d < ©k < 7r/2 

to arrive at H RF = (Vh 2 /U) + £ k ,± E u a l,± a K±- The 
eigenenergies of the ordered-state quasiparticles are given 
(for h > 0) by 
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The effective magnetic field which helps split the quasi- 
particle bands is fixed by solving the self-consistency 
equation h = Usq. For this simple model an explicit 
expression can be given for the right hand side and we 
obtain 



f(E^)-f(E+) 



= 1, 



(2.4) 



where f(E) is the zero-temperature Fermi function with 
the chemical potential determined by (1/V) Ylkifi^k ) + 



B. Specific toy model 

The above equations are valid for an arbitrary band 
dispersion relation ek. In Ref. free fermions with 



parabolic dispersion were considered. For this case states 
with spiral magnetic order always have higher energy 
than uniform magnetization (ferromagnetic) states. In 
the present paper, however, we want to construct a model 
in which the minimum (mean-field theory) energy oc- 
curs for a spiral state with wave vector Q* ^ 0, i.e., a 
model for which spiral magnetic order is favored over 
ferromagnetism. 16 This circumstance is achieved most 
simply by choosing a quasi-one dimensional model, i.e., 
by choosing 

W 

e k = — [1 - cos(k z a)} (2.5) 

for — ir/a < k z < 7r I a independent of k x and k v , where W 
defines the bandwidth. It is also convenient to consider 
the case of a half- filled band, i.e. a = n c ly/3 , where 
n c is the electron density. We study this pedagogical 
toy model in order to illustrate the general relationship 
discussed earliei^ between spiral magnetic order and spin 
supercurrents, not with the objective of modeling any 
specific material. In practice helimagnetism occurs for 
various reasons in several different types of materials; see 
for example Ref. and work cited therein. We will 
briefly discuss some of these materials in Sections IV and 
V. 

Choosing a quantization axis perpendicular to the 
spiral-order plane, the spin-projected current density is 
given by 

j' = ^E^< c L<W>. (2-6) 

The number operator average in the mean-field-theory 
state can be expressed in terms of Fermi occupation fac- 
tors for the quasiparticles of the Hartree-Fock Hamilto- 
nian leading to 

It is then straightforward to show that ek = £-k implies 
that j X = -j T . 

An alternative expression follows from Eq. (|1.2|l . which 
was discussed earliest and is explained in more detail in 
Appendix ^ The mean-field theory expression for the 
energy of a state with spiral-wavevector Q is 

<Q) 4 + /( ^ } + E * f(E ^ ■ (2 ' 8) 

k 

For all the numerical results discussed in this paper we 
explicitly checked the equivalence of Eqs. I|1.2(l and l|2.7|l . 
Note that a spin current can be present even though the 
quasiparticle population is in equilibrium. Elastic scat- 
tering from occupied to unoccupied quasiparticle states 
cannot provide a spin-current decay mechanism. The 
spin currents are instead carried collectively, as in the 
case of dissipationless charge transport in BCS super- 
conductors, and can decay only by collective processes!^ 
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FIG. 2: Quasiparticle bands E~£ (a) and (d) for two 
different values of Q. The interaction strength U is charac- 
terized by 2n c U/(jvW) = 1. For comparison, we also show 
the dispersion relations tk,+Q/2 (b) and tkz-Q/2 ( c ) for zero 
order parameter as dashed lines. 

that allow the phase, which specifies the spiraling mag- 
netization orientation, to slip. 

C. Numerical results 

Our numerical results are summarized in Figs. [2] and [3] 
We choose Q — Q z to be in the z-direction. In Fig. [21 we 
plot the quasiparticle energies E^ for Q = 0.7ir/a and 
Q = ir/a. 

The dashed curves represent £k,±Q/2 f° r spin-up and 
spin-down electrons, respectively; note that the bands 
are degenerate for k z a — n when the order parameter 
vanishes. The quasiparticle bands (solid lines) are 

split by ^(e k+Q / 2 - ek-Q/2) 2 + 4/i 2 , weakening the dis- 
persion of both the occupied band and the empty band. 
The self-consistently determined values of the order pa- 
rameter sq, the magnetic condensation-energy density 
e C ond = £o — e(<9) (where e = u e W(l/2 — 1/n) is the 
energy density of the state with h = 0), and the spin 
supercurrent density j = jj = — j± = are plotted as a 
function of the spiral wave vector Q in Fig. 01 For the 
half-filled band case we consider the condensation energy 
is maximized for Q = n/a. In accord with Eq. I|1.2|l no 
spin supercurrent is present at this value of Q. For other 
values of Q the derivative de(Q) / dQ is finite and the spin 
supercurrent is non-zero. In Fig. we have used dashed 
lines in the regions where dj(Q)/dQ is negative, to em- 
phasize that the spiral state is not metastable. Energy 
can always be gained in this regime without changing the 
total number of phase windings by introducing domains 
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FIG. 3: The order parameter sq normalized to the electron 
density n c , the magnetic condensation-energy density e con d 
normalized to the energy density eo of the disordered state, 
and the spin supercurrent density j — ji = —j± normalized 
to jo = en c Vk F as a function of the ordering wave vector Q 
for 2nJJ j '(ttW) = 1. The dashed lines indicate an instability 
regime against phase separation into regions with larger and 
smaller Q. 



with different phase winding rates. 

This example clearly shows that spiral magnetization 
configurations do not necessarily lead to spin supercur- 
rents; instead they occur when the system is driven from 
its ground state into a state with a non-optimal spiral 
wavevector Q. In the case of a fcrromagnet, i.e., a sys- 
tem for which the energy is minimized by Q = 0, it 
might* be possible to do this by driving the system from 
equilibrium using spin-selective transport currents sup- 
plied by ferromagnetic electrodes. This possibility has 
not yet been demonstrated experimentally, however, and 
may require combinations of material characteristics and 
geometries that are difficult to achieve. This difficulty 
motivates trying to find other strategies for realizing per- 
sistent spin currents. As we discuss in Section IV and 
V of this paper, in the case of a system with a spiral 
ground state persistent spin currents can be generated 
simply by applying an external magnetic field in the x- 
y plane. The in-plane field drives the spiral state to a 
soliton lattice state which does carry a persistent spin 



current. Before turning to this central portion of our pa- 
per, we first briefly comment on the relationship between 
our discussion of persistent currents in ferromagnets and 
Anderson's discussion^ of superconductivity in terms of 
magnetic order in an effective spin model. 



III. ANALOGY WITH BCS THEORY 

To understand the physics behind our results for the 
dissipationless spin currents it is instructive to invert 
Anderson's analogy and address magnetic properties in 
terms of superconductivity. An analogy between Joseph- 
son junctions and tunnel junctions between ferromag- 
netic metals has been presented recently22, which is in 
the same spirit as the present discussion. The spin super- 
currents discussed in this paper and the electrical super- 
currents supported by a Cooper pair condensate in BCS 
superconductors appear in a similar way when a particle- 
hole transformation is performed in one of the spin sub- 
spaces, say for spin down, to convert particle-particle 
order into particle-hole order. For example, the defini- 
tion of the order parameter in the microscopic model, 

• S Q = ( 1 /^)Ek( c k+Q/2,T Ck -Q/2a)' is mapped to the or- 
der parameter of a superconductor with a momentum Q 
pair condensate by this transformation. A finite value 
of the order parameter sq corresponds to the formation 
of "Cooper pairs" , consisting of an electron and a hole. 
These pairs carry no electric charge so that no charge su- 
percurrent arises from the order. The formation of pairs 
with a finite momentum is partly analogous to the pro- 
cess leading to Fulde-Ferrell-Larkin-Ovchinnikov states^ 
in a superconductor. 

We can achieve a qualitative understanding of our nu- 
merical results for the persistent spin current by com- 
bining Eq. qi.2p. Eq. (|2.3f) and Eq. <|2.8[) . and ignoring 
the dependence of the order parameter magnitude on Q. 
Then the contribution to the current from each wavevec- 
tor k in Ea. l2.8l depends only on the velocities of the elec- 
tron and the hole in the "Cooper pair", vj. = de^/d(Kk) 
and v_k = — Vk- (These bare velocities can be identified 
with the slopes of the dashed lines in Fig. [21) At finite Q, 
however, the velocity of a "Cooper pair" is finite and the 
sum over all wavevectors k in general produces a result 
that is not zero. Thinking in this way we can see directly 
that the point Q = ir /a, for which we found in the pre- 
vious section that the condensation energy is maximal 
and the total spin supercurrent is zero, is special. The 
reason is that the contribution of each "Cooper pair" is 
neutralized by that of another pair with reversed wave 
vectors for both electron and the hole. 



IV. PERSISTENT SPIN CURRENTS IN 
HELIMAGNETS 

Two classes of helimagnets have been extensively stud- 
ied in the past. The first is comprised of MnSi and related 
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materials that have long-period helical ground states 
because of the Dzyaloshinskii 22 instability of ferromag- 
netism in systems without inversion symmetry. These 
systems have been studied extensively very recently 2 ^ be- 
cause they provide an example of an itinerant-electron 
magnetic system in which the ordering temperature can 
conveniently be driven to zero by applying pressure, and 
because they appear to show non-Fermi-liquid behavior 
associated with this nearby quantum critical point. They 
also have the potential advantage for the phenomenon 
of interest here that the chirality degeneracy, present in 
the microscopic model studied in Section II for example, 
is lifted— by broken inversion symmetry. In this case 
however, the plane of the spiral magnetic order is fixed 
only indirectly by anisotropy in the gradient term in the 
Landau-Ginzburg energy functional and the easy plane 
anisotropy that we require is rather weak. 

The second class of materials consists of the heavy met- 
als Tb, Dy, and Ho which have helical ground states due 
to frustration induced by RKKY interactions between 
the rare-earth moments and do have strong easy-plane 
anisotropyiSSiSSiSl At present, it appears to us that these 
are the most promising materials for the realization of 
persistent spin currents. 

In the case of a superconductor, charge supercurrents 
can be generated quite simply by biasing the sample with 
an external current source. In the case of an easy-plane 
ferromagnet, spin supercurrents can be generated by bi- 
asing the systems with an external spin current. One 
possible scheme for realizing such a bias using four fer- 
romagnetic contacts is discussed in Ref. |!£ but its suc- 
cess depends on avoiding spin-flip processes at the inter- 
faces and on having thin-film samples with spin-diffusion 
lengths that are larger than the sample size. For systems 
with a spirally-ordered ground state, however, there is an 
easier and more straightforward way of generating per- 
sistent spin currents which does not depend on achieving 
spin-current biasing. Instead, persistent spin currents 
can be generated in the ground state of the system sim- 
ply by applying a magnetic field in the plane defined by 
the spiral order parameter. In the continuum model that 
we study, the persistent-current state in a helimagnet is a 
soliton-lattice state for which the lattice Bloch wavevec- 
tor is not a good quantum number. For this reason we 
cannot easily describe these states microscopically. We 
therefore use a phenomenological model for the following 
discussion. 



A. Magnetic order in the presence of an in-plane 
magnetic field 

A magnetic field in the x-y plane alters the spiral 
state. At sufficiently large field strengths, it is clear 
that the magnetic field will lead to nearly uniform spin- 
polarization along the field direction. As we show, for 
fields below a critical strength, the spiral state is driven 
by the field to a soliton-lattice state, in which the mag- 



netization orientation slips periodically from the field- 
direction orientation. 

We assume that due to easy-plane anisotropy the mag- 
netization is in the x-y plane and has a spatially con- 
stant amplitude. The orientation variation of the mag- 
netization density m(r) = gp^ ( s ( r )} is parametrized by 
the angular variable d(r) via m(r) = mo[cos(^(r))x — 
sin(-i?(r))y]. For clarity, we choose Q* in z-direction, 
Q* = (0, 0, Q*) so that only depends on the z-position. 
However, the following discussion does not depend on 
this specific choice. To determine the magnetic order at 
T = we have to minimize the energy density 



dz 

T~ z 



■\d z d{z)-Q*\ 2 - B- m(z) 



(4.1) 



where L z is the length of the sample in z-direction, the 
spin stiffness p s characterizes the energy cost for spiral 
order with wavectors close to the minimal value Q*, and 
B is the in-plane magnetic field. 

Without loss of generality we choose the inplane mag- 
netic field along the x direction, B = (B,0, 0). For this 
particular choice of directions, $ is the angle between the 
local magnetization direction and the external magnetic 
field. This model for a spiral state in an external field 
(Eq. I|4.1|)l is equivalent to the Pokrovsky-Talapov (PT) 
model, reviewed for example by P. Bak(2£ used originally 
to model commensurate-incommensurate transitions and 
more recently to model the influence of an in-plane field 
on spontaneous coherence broken symmetry states in bi- 
layer quantum Hall systems. In using this model we as- 
sume that the energy density depends on the gradient 
of the orientation angle ?? and that, for qualitative pur- 
poses, we can expand around the value of the gradient 
that minimizes the energy density. This chiral model fa- 
vors a particular sign for c? z $ and does not account for 
the symmetry-based expectation that the energy density 
will have minima for d z d = ±Q* . Hence it accounts only 
for the effect of a magnetic field on a spiral state domain 
with a particular chirality; we speculate later on persis- 
tent spin currents near boundaries between domains with 
opposite chirality. A more general model of the gradient 
energy that does not favor a particular chirality would 
have the form f(dz/L z )(p s /2)[-(d z tf) 2 /2+(d z tf) 4 /<iQ* 2 } 
plus a constant. Finally, our continuum model ignores 
the possibility of locking to configurations that are com- 
mensurate with the underlying lattice. Locking to com- 
mensurate configurations has been studied for helimag- 
nets in fields for both non-chiral2£ and chiral 3 ^ models. 
Indeed, magnetoelastic coupling is known to play an im- 
portant and complex role-i in the response of helimag- 
nets to external magnetic fields. Our use of a chiral con- 
tinuum model for the following discussion certainly does 
not capture the full richness of real helimagnets in exter- 
nal magnetic fields. Our intention here is to demonstrate 
simply the general property that external fields induce 
persistent spin currents within each chiral domain, using 
a model that is well understood 28 31 . In the following 
discussion we use well known properties of the classical 
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FIG. 4: Energy versus in-plane magnetic field for our spiral 
state model. The soliton-lattice state minimizes the model 
energy for B < B c , while the state in which the magnetization 
is aligned with field minimizes the energy for B > B c . The 
soliton lattice state energy approaches the spiral state energy 
as B — > 0. The energies are normalized to B c mo. 



sine-Gordon model to evaluate the persistent spin current 
within a chiral domain of a helimagnet. 

Minimization of the energy with respect to variations 
in $(z) leads to the sine-Gordon equation 



d 2 i){z) 
dz 2 



2 sm(v)(z)) 



(4.2) 



where £ = J p s /Bmo- For fields B larger than a crit- 
ical value B c , the ground state magnetization is uni- 
form with = 0. For B just below B c the energy 
can be lowered by incorporating isolated 2ir solitons; 
$ss{z) = 4 arctan[exp(±(z — zo)/£)] for a soliton centered 
at zq. At B c , the uniform state has the same energy as a 
state which accommodates a single soliton. By compar- 
ing energies we find that in terms of the notation defined 
above 



16 mo 



(4.3) 



For B < B c the magnetic state can be described as a 
soliton lattice (SL) with period a. The associated wave 
vector Q = 2-k ja varies from Q = at B = B c to Q = Q* 
at B = 0. In the latter limit the soliton-lattice state ap- 
proaches the spiral magnetic order state, discussed from 
a microscopic point of view in the first part of the pa- 
per. The total phase change along the sample is given by 
QL Z . 

The results collected in the following paragraph are 
adapted from Ref.|j^]and earlier work cited therein. Min- 
imization of the model's energy in the soliton lattice state 
gives, cos(i?(z)/2) = — sn[(z — zo)/r)£,ij\ where sn denotes 
the sine-amplitude Jacobian elliptic function^ and rj is a 
constant which depends on the strength of the magnetic 
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(4.4) 



where E(rj) = J^ 2 dj3\J\ — rj 2 sin 2 (3 denotes the com- 
plete elliptic integral of the second kind. 32 From that we 
see that the limit rj — 1 applies at the SL state with 
B = B c , and that rj = describes the opposite situation 
of a vanishing in-plane magnetic field B = 0. The wave 
vector Q of the SL is given by 



Q_ (tt/2) 2 
Q* K( V )E(rj) ' 



(4.5) 



where K (77) = J d(3/ y 1 — rj 2 sin (3 denotes the com- 
plete elliptic integral of the first kindi 3 ^ We note that 
Q = for B = B c , and Q = Q* for 5 = 0. Finally, 
the expression for the energy density of the SL, which we 
relate to persistent spin currents below, is 



(4.6) 



The SL energy density is plotted as a dashed line in 
Fig. J3J. For B < B c the SL is energetically favored 
over a uniform state, whose energy is (l/2)p s Q* 2 — Brno. 
For B ^ the SL lattice is always favored over the spiral 
state (with wave vector Q*). The energy of the latter 
state is the zero of energy of our model. 

A cautionary note is appropriate at this point. We 
do not expect that the expansion of the gradient energy 
around the spiral wavevector that appears in our model 
will be realistic for typical spiral magnets all the way 
down to the Q = state that is favored by the in-plane 
field. It follows that our model is most realistic physically 
in the B < B c region. 



B. Spin supercurrent 

We now demonstrate how the in-plane magnetic field 
gives rise to a persistent spin current. To derive an 
expression for the spin current density we introduce a 
spin-dependent, spatially constant vector potential A CT = 
A a i. The energy density of the system then reads 



1 v ^ f dz 
2^./I7 



2 [d- Z - Q +a Yc A 



-Bm cos $ 



he 



(4.7) 



Note that here we have explicitly used the coincidence of 
r) with the angle between the magnetization direction and 
the magnetic field; for a different orientation of the mag- 
netic field in the x-y plane, Eq. I|4.7|l adopts a different 
form. The spin-dependent current density follows from 
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FIG. 5: Spin supercurrent density j = jf induced by an in- 
plane magnetic field. In the spirally ordered state at B = 
there is no current and for B > B c the current is constant. 



FIG. 6: Sketch of the setup discussed in the text. The direc- 
tion of the magnetization is indicated by the arrows. In the 
wires it is parallel to the z-axes along the transport direction, 
while in the helimagnet (shaded grey) it is circulating in the 
x-y plane as shown in the lower part of the sketch. 



the derivative = cde(Af, A±)/dAf at Af = A± = 0. 
Making use of the relation J dz sm{d(z)) = 0, that holds 
for all soliton lattices with arbitrary wave vector, we find 
that 

jt = -ft (Q - Q*) = -h ■ ( 4 -8) 

At B = we have a spiral state with Q = Q* , and the 
spin current vanishes. With increasing strength of the B 
field the wave vector Q decreases and a finite spin super- 
current arises. At B = B c the spin current jj reaches 
its maximal value jo = ~{e/fi)psQ* ■ The limit of small 
magnetic fields, B <C B c , is described by 

The spin current density as a function of the in-plane 
magnetic field is shown in Fig. [3J 

In conclusion, the SL state of a helimagnet carries a 
finite persistent spin current, without any electrical cur- 
rent, when a magnetic field is applied in the plane defined 
by the spiral order. 



V. DISCUSSION 

We have shown in this paper that persistent spin 
currents can be induced in helimagnets by applying a 
magnetic field in the plane defined by the spiral order. 
It is interesting to speculate on the possibility of ex- 
ploiting this effect for new types of spintronic devices. 
Current spintronic devices make use of magnetotrans- 
port effects like anisotropic magnetoresistance, giant 
magnetoresistance?^ and tunnel magnetoresistance^i all 
of which follow from the dependence of quasiparticle 
transport on the collective magnetic state. The effect 
we have discussed here is distinctly different in that it 



is an equilibrium spin current that is carried collectively 
rather than by quasiparticles. 

To illustrate the possibility of realizing interesting 
magnetotransport effects based on field-induced persis- 
tent currents in helimagnets, consider the case illustrated 
in Fig. |H| We consider current flow along a ferromagnetic 
wire which contains a helimagnet. The magnetization 
direction in the ferromagnets is assumed to be aligned 
along the direction of the wire by magnetostatic energy, 
and the electrical current in the ferromagnet will be spin- 
polarized along the direction of the magnetization. We 
also assume that the helimagnet is grown so that it has its 
wavevector along the wire direction and the spiral mag- 
netization is in the perpendicular plane; this would be 
the case for example for rare earth metal single-crystal 
helimagnets with their hep c-axis aligned with the wire. 
The quasiparticles of the helimagnet would then tend to 
have their spins aligned in the plane perpendicular to 
the wire axis, suppressing their ability to carry currents 
that are spin-polarized along the wire axis, and increas- 
ing the resistance of the overall system. In the presence 
of a field perpendicular to the wire axis, however, the 
helimagnet would have a persistent spin current polar- 
ized along the wire axis, eliminating the need to carry 
the spin current with quasiparticles. We predict that be- 
cause of the persistent spin currents, an in-plane field 
could therefore alter the resistance of this system. Since 
the persistent spin current has a sign that depends on 
the chirality of the spiral, the sign of the resistance de- 
pends on this quantity. The effect should be strongest in 
samples with a single chiral domain. 

This same effect could also be of interest in a simi- 
lar setup where the helimagnet is connected to normal 
metallic leads. Then the spin current in the helimag- 
net results in a different resistance for the electron spins 
entering from the leads. In this fashion the helimagnet 
could act as a tunable spin filter which is controlled by 
an external magnetic field. 

It may also be possible to realize spin torque^ effects 
in helimagnets that are distinct from those which occur 
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in ferromagnets. Spin torque occurs when the spin po- 
larization of a quasiparticle current changes with posi- 
tion by altering the spin distribution of current-carrying 
states, transferring the spin carried by the quasiparticle 
to the collective magnetic coordinate. The field-induced 
persistent spin current in a helimagnet will have a strong 
spatial dependence, a sign change in fact, near a bound- 
ary surface between domains with opposite chirality. The 
divergence of this spin current makes a contribution to 
the time-dependent spin polarization that is cancelled by 
the torque due to the external magnetic field when the 
domain wall is in equilibrium. When the domain wall is 
not in equilibrium, for example by being displaced from 
a pinning center, these two contributions will not cancel 
and the net torque will lead to a time-dependent magne- 
tization orientation. Conversely, the equilibrium position 
of a pinned domain wall will be sensitive to the presence 
of persistent currents and will therefore be altered by an 
in-plane field. 

In conclusion, we have studied the relation between spi- 
ral magnetic order and dissipationless spin transport. We 
have demonstrated the possibility of dissipationless spin 
currents in states with equilibrium quasiparticle popu- 
lations. These spin currents are always associated with 
spiral magnetic order. On the other hand, the existence 
of spiral magnetic order is, in general, not sufficient to 
guarantee collective spin transport. In particular heli- 
magnets, which have spiral ground states, do not sup- 
port persistent spin currents. However, persistent spin 
currents can be induced in these systems simply by ap- 
plying an external magnetic field. 
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APPENDIX A: PROOF OF EQ. (2) 

In order to prove that the spin-current density is 
related to the derivative of the total energy with re- 
spect to the wave vector Q as given in Eq. (|I.2() we 
extend the Hartree-Fock Hamiltonian Eq. (|2.1|) by for- 
mally introducing an additional parameter Q via the 
replacement e k+Q/2 -> e k+( Q + Q)/ 2 and e k _Q/ 2 -> 
e k-(Q+Q)/?- O ur original Hamiltonian Eq. (|2.1|) corre- 
sponds to H HF (Q, Q)|q =0 - The quasiparticle ground 

state for given Q and Q is an eigenstate and, therefore, 
the Hellmann-Feynman theorem applies, 



d{H (Q, Q)) 



<9Q 



dH Mt (Q, Q) 



(Al) 

Q=0 

we find that 



Q=o \ <9Q 
By comparing the r.h.s with Eq. I|2.6[) 
Eq. (fAlfl e quals to V(7i/2e)(jf - jj,). To evaluate the 
l.h.s of Eq. iJMJ we observe that the quasiparticle ener- 
gies for H MF (Q, Q) as well as the value of h depend 
on Q and Q only in the combination Q + Q. As a con- 
sequence, we can replace d/dQ by d/dQ. Together with 
jj = — j;, this immediately proves Eq. JT2J. 
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